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The effect of electronic interactions in graphene with vacancies or resonant scatterers is inves- 
tigated. We apply dynamical mean- field theory in combination with quantum Monte Carlo sim- 
ulations, which allow us to treat non-perturbatively quantum fluctuations beyond Hartree-Fock 
approximations. The interactions narrow the width of the resonance and induce a Curie magnetic 
susceptibility, signaling the formation of local moments. The absence of saturation of the suscepti- 
bility at low temperatures suggests that the coupling between the local moment and the conduction 
electrons is ferromagnetic. 



Introduction. Since its isolation^^^, single layer 
graphene has attracted a great deal of attention, due to 
novel features. The massless nature of the charge carriers 
implies that the density of states vanishes at the Fermi 
energy in a neutral layer—. The low density of states in 
its vicinity allows for the formation of sharp resonances 
due to vacancies^ or impurities like hydrogen which form 
a strong covalent bond with the carbon atoms, the so 
called resonant impurities^. These resonances have been 
observed in graphite surfaces^. 

The enhancement in the density of states by those res- 
onances and the electron electron interaction favor the 
formation of local moments. The states associated to the 
resonances differ from those induced by coupled magnetic 
dopants in a number of ways: i) They are built up from 
the same tt orbitals as the conduction band of graphene, 
ii) The resonance state is orthogonal to the conduction 
states, and the hopping between the resonance and the 
extended states vanishes, and iii) They extend over a 
large region near the defect, as there is no gap in the 
spectrum to confine them. 

Mean field arguments based on the enhancement of 
the local density of states favor the formation of a 
static magnetic moment, as shown in a number of 
calculations^"— . It is known that quantum fluctuations, 
not included in Hartree-Fock approximations, determine 
the low temperature properties of magnetic impurities 
in solids, described by the ferro- or antiferromagnetic 
Kondo model^^. The differences between the resonant 
levels in graphene and magnetic impurities in metals and 
graphenei^ — imply that mean field results cannot be ex- 
trapolated in a straightforward way to low temperatures. 

In the following, we analyze the electronic properties of 
graphene by non perturbative methods beyond a static 
mean field approximation. Our results show that, for 
reasonable interaction strengths, the main features are 
well described by assuming the existence of a fluctuat- 
ing magnetic moment. This moment is not quenched 
at the lowest accessible temperatures, suggesting a fer- 
romagnetic coupling with the conduction band. For a 
finite concentration of resonant impurities, the absence 
of competition between the antiferromagnetic Kondo ef- 
fect and the RKKY interaction might lead to ferromag- 




FIG. 1. (Color online). Sketch of the model studied in the 
text. The presence of a vacancy, or a resonant scatterer, pas- 
sivates one of the lattice sites. The resulting localized state 
near the Dirac energy becomes spin polarized. 



netism, provided that the concentration of impurities is 
large enougbii — . 

The model. We describe the tt band of graphene by 
a nearest neighbor tight binding model, with a hopping 
parameter t. The effect of a resonant scatterer strongly 
bound to a given site is taken into account by shifting 
the on-site energy by an amount eo- In the limit |eo| ^ t 
the model describes an unrelaxed vacancy. For |eo| ^ t, 
a resonance near the Dirac energy builds up. This reso- 
nance moves to the Dirac energy and its width vanishes 
in the limit |eo|/t oo. The main features of the model, 
including the possibility of a local moment at the reso- 
nance, are shown in Fig. [TJ 

We assume that the long range part of the electron- 
electron interaction is screened, and we describe the elec- 
tron electron interaction by an on-site Hubbard repulsive 
term, U . The full hamiltonian is 



{i,j)cr cr 



(1) 



where c[^J annihilates (creates) an electron with spin 

a =t, i on lattice site i and where riia = cl^c-^ denotes 
the corresponding number density. The hopping t is only 
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FIG. 2. (Color online). Sketch of the lattice system with two 
sub- lattices A and B (left panel). The impurity is denoted 
by an open circle. We map the lattice on a two-site cluster 
consisting of only one unit cell (right panel). The cluster is 
embedded in an interacting medium which determined by a 
local self-energy S. The self-energy is calculated either by 
DMFT or by second-order perturbation theory. 



finite for neighboring lattice sites (denoted by (i^j)). The 
lattice site i = corresponds to the impurity where the 
local on-site energy eo is non-zero. We study the model 
at half filling, i.e., at chemical potential /i = 0. 

We approximate the interacting lattice problem by as- 
suming that the impurity site where cq ^ and a neigh- 
boring site are attached to an effective medium described 
by a local self-energy S as sketched in Fig. [2l The effect 
of the impurity on the bath vanishes in the thermody- 
namic limit and can therefore be neglected^. The re- 
sulting two-site cluster problem is solved numerically by 
quantum Monte Carlo (QMC) simulations. We employ 
QMC methods in continuous imaginary time which per- 
form a systematic expansion in the interaction term of 
the Hamiltonian^i^. The QMC solution of the cluster 
problem is numerically exact and fully incorporates inter- 
actions and quantum fluctuations. The analysis can be 
regarded as the solution of a quantum impurity problem 
where the interaction effects are included at the impurity 
site and at a close neighbor, and approximated by means 
of a local self-energy in the surrounding medium. 

In order to calculate the self-energy necessary to de- 
termine the effective medium, we start by simulating the 
homogeneous lattice system with eo = 0. We calculate 
the self-energy necessary to obtain the effective medium 
using either Dynamical Mean Field Theory (DMFT)^^ or 
second-order perturbation theory in U /t. DMFT fully in- 
corporates quantum fluctuations local to the cluster but 
ignores spatial fluctuations. The medium depends on the 
self-energy of the cluster system and has to be calculated 
self-consistently by an iterative procedure using QMC 
simulations. Subsequently, the impurity is added to the 
system and one additional QMC simulation is performed 
using the medium of the converged DMFT calculation. 

In order to check the quality of the DMFT approxi- 
mation, we additionally calculate the self-energy of the 
homogeneous lattice system using second-order pertur- 
bation theory in U/t. This self-energy - instead of the 
self-consistently determined DMFT solution - is used to 
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FIG. 3. (Color online). Density of states D{uj) at the site next 
to the impurity at /3t = 10 calculated by DMFT in combina- 
tion with QMC. The analytic continuation was performed by 
maximum entropy. Top: density of states for e^/t — 3 and 
different values of U/t. Bottom: density of states for U/t = 4 
and different values of eo/t. The insets show a blowup of the 
region near the Dirac energy. 



calculate the bath of the impurity problem which is then 
again solved by QMC. The perturbative self-energy does 
not include quantum fluctuations to all orders. How- 
ever, it incorporates non-local effects of the actual lattice 
structure which are neglected by DMFT. Thus we are 
able to test the influence of non-local correlations and 
the accuracy of the DMFT approximation. 

QMC methods map quantum-mechanical systems on 
a classical one at the expense of an additional dimension 
which - in most cases - is an imaginary time dimension. 
Thus, QMC can only provide dynamical data for imag- 
inary times or frequencies. The necessary analytic con- 
tinuation to physically relevant real times or frequencies 
is usually performed by maximum entropy techniques'^. 
We use a standard maximum entropy implementation^^ 
to calculate the interacting density of states for real fre- 
quencies. 

Results. The density of states at the site next to the 
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FIG. 4. (Color online). Density of states at the site next to the 
impurity at /3t = 10 calculated by second- order perturbation 
theory. Parameters are /3t = 10, U/t = 1, and different values 
of eo/t. The inset is a blowup of the region near the Dirac 
energy. 



impurity is shown in Fig. [3l We find a resonance whose 
width decreases either by increasing eo/t or by increasing 
U/t. The reduction in width by the interactions is a 
characteristic feature of magnetic impurity problem^i^. 
It indicates the decoupling of the impurity degrees of 
freedom from the conduction band. For eo > 2t and 
U > 2twe find a second resonance to the left of the Dirac 
point which is shifted further to the left for increasing 
values of eo/t or U/t. There appear also satellite peaks 
at high energies. These peaks are related to transitions 
involving configurations where the resonance hosts zero 
or two electrons. The positions of these peaks are shifted 
by an amount proportional to eo, confirming that they 
are related to the impurity. 

We have repeated the previous calculations using an 
input self-energy obtained from second-order perturba- 
tion theory. The results are shown in Fig. |4] and are 
consistent with those shown in Fig. [3] for U/t < 3. One 
clearly sees the formation of a resonance that becomes 
sharper and shifts towards uj/t = for increasing eo/t. 
We also have a peak to the left of u/t = that is shifted 
with increasing eo- 

We calculate the magnetic susceptibility of the impu- 
rity in the imaginary-time framework of the QMC via 



Xs 



/5 



dr(K(r)-n;(r))(nt(0)-n;(0))) 



(2) 
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where r denotes imaginary time, ncr(r) = e 
and j3 = l/ksT. As usual, T denotes temperature and 
kB Boltzmann's constant. Fig. [5] shows xs as function 
of the inverse temperature. A fully localized spin pos- 
sesses the susceptibility xioc = Pt. The numerical results 
are consistent with a Curie dependence on temperature, 
X oc 1/T, suggesting the formation of a local moment. 
This moment is not fully localized at the positions clos- 
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FIG. 5. (Color online). Local magnetic susceptibility xs as 
function of the inverse temperature calculated by DMFT in 
combination with QMC. Top: susceptibilities for e^/t — 3 and 
different values of U/t. Bottom: susceptibilities for U/t = 4 
and different values of eo/t. 



est to the impurity, and the measured xs is smaller than 
~ Xioc/2. As the interaction increases the magnetic mo- 
ment becomes better defined and more localized near the 
impurity. We find no sign of saturation of the suscepti- 
bility down to the lowest temperatures. This result is 
consistent with a ferromagnetic Kondo coupling. 

Discussion. We have studied the effects of interactions 
in the presence of a resonance in the graphene electronic 
spectrum. We use a local approach, and fully include 
quantum fluctuations. The interaction is described by a 
local Hubbard term, and we do not consider the imperfect 
screening in graphene near the neutrality point which is 
expected in suspended layers. Given this interaction, our 
calculation can be regarded as the solution of a quantum 
cluster problem. The effect of interactions in the medium 
surrounding to the cluster is described by means of a local 
self-energy. The results obtained when this self-energy 
is calculated by DMFT and by perturbation theory are 
mutually consistent, in the regime where perturbation 
theory is valid. The off diagonal corrections to the self- 
energy not included here are mostly due to the long range 
part of the interaction^i. They lead to an increase in the 
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Fermi velocit}*^^. The density of states at low energies is 
reduced, favoring further the formation of local moments. 

We find that interactions reduce the width of the res- 
onance induced by the impurity, and lead to a mag- 
netic susceptibility which grows at low temperatures as 
Xs oc 1/T. The local moment is not quenched at the low- 
est temperatures studied. This is consistent with the ex- 
istence of a ferromagnetic coupling between the moment 
and the valence electrons. Note that the exchange mech- 
anism which leads to an antiferromagnetic interaction for 
magnetic impurities coupled to a conduction band does 
not exist in the case of a resonance built up from the or- 
bitals which also give rise to the conduction band. The 
existence of the vacancy does not lead to virtual hoppings 
between extended and localized states. An electron oc- 
cupying the resonant state interacts with a conduction 
electron through the onsite Hubbard term. This cou- 
pling favors a ferromagnetic alignment of the spin of the 
electron in the resonance and the spin of the conduction 
electrons. 

The scale at which the local moment studied here leads 
to significant effects depends on the value of which is 



not very precisely determined in graphene. Calculations 
based on the Local Density Functional Approximation 
suggest^ U/t ~ 1, while quantum many body calcula- 
tions for aromatic molecules give^ — U/t ~ 3. The value 
of U/t is bounded by U/t ~ 4.5, above which graphene 
should become antiferromagnetic^. 
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